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channel driven by longitudinal flow 

Kaifu LucE 

CAS Key Laboratory of Soft Matter Chemistry, Department of Polymer Science and Engineering, 
University of Science and Technology of China, Hefei, Anhui Province 230026, P. R. China 

Ralf MetzleiQ 

Physics Department, Technical University of Munich, 85H1 Garching, Cermany and 
Physics Department, Tampere University of Technology, FT33101 Tampere, Finland 

(Dated: January 15, 2013) 

Using analytical techniques and Langevin dynamics simulations, we investigate the dynamics 
of polymer translocation into a narrow channel of width R embedded in two dimensions, driven 
by a force proportional to the number of monomers in the channel. Such a setup mimics typical 
experimental situations in nano/micro-fluidics. During the the translocation process if the monomers 
in the channel can sufficiently quickly assume steady state motion, we observe the scaling r ~ N/F 
of the translocation time r with the driving force F per bead and the number A'' of monomers per 
chain. With smaller channel width R, steady state motion cannot be achieved, effecting a non- 
universal dependence of r on A'^ and F. From the simulations we also deduce the waiting time 
distributions under various conditions for the single segment passage through the channel entrance. 
For different chain lengths but the same driving force, the curves of the waiting time as a function 
of the translocation coordinate s feature a maximum located at identical Smax, while with increasing 
the driving force or the channel width the value of Smax decreases. 



I. INTRODUCTION 

The transit of biopolymeric chains across membranes 
through nanopores, the translocation process, is a recur- 
rent theme in ceh biology [H-Q . In biological membranes 
the nanoporc is typically constituted by a channel pro- 
tein, the transmembrane proteins of the haemolysin fam- 
ily being a prominent example (Bf. One family member, 
a-hacmolysin, is also often employed in single nanoporc 
setups in vitro @ . More recently, solid state nanopores in 
artificial supports such as silicon compound membranes 
are manufactured by ion or electron beam techniques and 
open up the possibility for controlled technological appli- 
cations 

In biology, important examples for translocation pro- 
cesses are the passage of RNA or proteins through pores 
in the nuclear membrane, of proteins across mitochon- 
drial and chloroplast membranes or through the endo- 
plasmatic reticulum, the cell-to-cell exchange of DNA 
across the cell walls, or the viral injection of RNA and 
DNA molecules [i| . In technology, biopolymcr transloca- 
tion is envisaged to be useful for rapid DNA sequencing, 
gene therapy, and, ultimately, towards controlled drug 
delivery. 

The transport of biopolymcrs through a nanoporc has 
attracted broad interest in the statistical physics commu- 
nity, as it represents a challenging problem in polymer 
physics [9I-I22I [2^ [25I . A quantity of particular inter- 
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est is the average translocation time r as a function of 
the chain length TV, usually assumed to follow a scaling 
law r ~ iV". The scaling exponent a hereby reflects the 
efficiency of the translocation process. Thus, for com- 
pletely directed (ratcheted) motion, one would expect 
T TV, while for normally diffusive unbiased transloca- 
tion T ^ N"^ 0, [1^ . Generally, the value of a differs from 
these limiting behaviors, as effected by the entropic de- 
grees of freedom of the chain-to-be-translocated. In par- 
ticular, when the translocation dynamics is governed by 
unbiased anomalous diffusion [26[ of the form (s^(t)) ^ t^ 
(0 < /3 < 1) in terms of the translocation coordinate s, 
the value of the scaling exponent a = 2//3 becomes larger 
than 2 0,113. 

During the passage of a long flexible chain through 
a nanopore, its two extremities cannot access the vol- 
ume opposite the membrane, and the monomer(s) in the 
pore are more of less unable to move. This reduction of 
the accessible degrees of freedom of the polymer involves 
a considerable entropic barrier during the pore passage. 
While unbiased translocation has been argued to feature 
the same scaling exponent a ^ l-\-2v ss, the free diffusion 
of the polymer by its radius of gyration [T2| - [l^ . where 
the Flo ry exp onent \s v = 0.588 in 3D and V2d = 0.75 
in 2D [23-[30|, the associated prefactor is much larger 
p^ . Thus efficient polymer translocation requires the 
presence of driving forces. Most experimental [Q-Sl and 
theoretical [H-Ill [13, [H El] studies focus on the driv- 
ing force provided by an external applied electric field, 
which mainly falls off across the pore. In this case, a de- 
pends on the translocation dynamics. For slow translo- 
cation, i.e., under low driving force and/or high friction, 
a « l + v « 1.588 in 3D 0,[3al, while a w 1.37 in 3D [Hi 
for fast translocation due to the highly deformed chain 



conformation on the trans side, reflecting a pronounced 
non-equilibrium situation. In addition, for transfoca- 
tion achieved by pulHng one end of the polymer with 
a constant force, a ~ 2.0 was established [Mill. The 
above examples of translocation involve constant forces 
throughout the process. A different situation is encoun- 
tered for translocating chains whose passage is effected 
by binding proteins, so-called chaperones or chaperonins 
3J) . Due to the fact that the chaperones bind more than 
one monomer of the chain-to-be-translocated, interesting 
variations occur in the effective driving force due to the 
"parking lot effect" . 

Here we consider a case motivated by state-of-the-art 
fluidic nanochannel setups, in which the driving force 
grows constantly with the progress of the translocation 
process, until a saturation is reached when the entire 
chain is in the channel. In such setups one or multiple 
nanochanncls branch off a micron-sized feeder channel 
(access hole) [35, . i36| . The driving force is provided by 
the liquid flow: while the flow fleld is quite low in the 
feeder channel, the flow velocity is much higher inside 
the narrow channel. Once one end of the chain enters 
the nanochannel, it is veritably sucked inside, the overall 
drag force on the chain increasing as more and more of 
the monomers enter the channel. Similarly, one might 
think of setups in which an electric fleld falls off in a long 
channel. If each monomer of the chain-to-be-translocated 
carries a net charge, the force acting on the chain will in- 
crease with the proportion of the chain inside the chan- 
nel. This is not the case for the relatively short channel 
[37| . where the driving force cannot be proportional to 
the translocated monomers once the flrst monomer has 
passed through the pore. In many cases such fluidic se- 
tups are built in a pseudo two-dimensional geometry, i.e., 
the channel is sandwiched between two parallel walls, 
that are close to each other. Having such scenarios in 
mind we consider the geometry depicted in Fig. [TJ We 
investigate this problem using analytical techniques and 
Langevin dynamics (LD) simulations. The paper is orga- 
nized as follows. In Section II, we describe our model and 
the simulations technique. In Section HI, we present and 
discuss our results. Finally, the Conclusions are drawn 
in Section IV. 



II. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard- Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential [13 • Excluded volume interactions be- 
tween beads are taken into consideration by a short range 
repulsive LJ potential 



(1) 




FIG. 1: Schematic representation of polymer translocation 
into a long 2D channel of width R. The driving force acting 
on the bead in the channel can be provided by an applied 
intra-channel electric field or a uniform flow of the solvent in 
the channel. We assume that the net force on the chain is 
proportional to the number of monomers, s{t), in the chain. 



between neighboring beads is modeled as a FENE spring 
with 



FENE 



(r) 



__fci?^ln(l-_ 



(2) 



where r is the distance between consecutive beads, k is 
the spring constant, and Rq is the maximum allowed sep- 
aration between connected beads. 

We consider the geometry shown in Fig. [1] two strips 
with separation R consisting of stationary particles with 
a bead-bead distance a form the walls of the channel, as 
well as the "membrane" containing the pore. Between 
all bead-wall particle pairs, there exists the same short 
range repulsive LJ interaction as described above. In the 
LD simulations, each bead is subjected to conservative, 
frictional, and random forces, respectively, with [s^ 



-V(C/lj + t/pENE) 



(3) 



for r < 21/6(7 and for r > Here, a is the diameter 

of a bead and e is the potential depth. The connectivity 



where m is the bead mass, ^ is the friction coefflcient for 
a single bead, v.; = is the bead velocity, and Ff is the 
random force satisfying the fluctuation-dissipation theo- 
rem. The external force is expressed as Fext — Fx, where 
F is the external force strength per bead, exerted on the 
translocating beads located inside the channel, and i is a 
unit vector in the direction along the channel. Therefore, 
the driving force for the whole chain is proportional to 
the number of translocated beads. Experimentally, this 
external driving force acting on the beads in the channel 
can be provided by an applied intra-channel electric fleld 
or a uniform flow of the solvent in the channel. 

In the present work, we use the LJ parameters s and 
tr and the bead mass m to fix the energy, length and 
mass scales respectively. The time scale is then given 
by tL.i = (to(t^/£)1/^. The dimensionless parameters in 
our simulations are i?o = 2, fc = 7, ^ = 0.7 and ksT = 
1.2 unless otherwise stated. The Langevin equation is 
integrated in time by a method described by Ermak and 
Buckholz lii in 2D. 
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Initially, the first bead of the chain is placed just in- 
side the channel (at x = 0.75, y = 0), while the remain- 
ing beads are under thermal collisions described by the 
Langevin thermostat to obtain an equilibrium configu- 
ration. The translocation time is defined as the time 
duration between the beginning of the translocation and 
the last monomer entering into the channel. Note that, 
in contrast to many translocation setups, we allow the 
chain to escape from the channel back towards the cis 
side of the membrane. Typically, we average our data 
over 1000 independent runs. 

Here, we should mention that the effect of hydrody- 
namic interactions is neglected in our simulation. Recent 
lattice Boltzmann [4l|, and molecular dynamics [i^ 
simulation results show that hydrodynamics is screened 
out in a narrow pore, which is the case here and in the 
experiments. 

III. RESULTS AND DISCUSSION 

In this Section, we present our simulation results and 
discuss them in view of scaling arguments provided in the 
following Subsection. 

A. Scaling arguments 

In standard translocation models, where the driving 
force acts solely in a short channel, only a single bead 
currently inside the channel experiences the driving force. 
For sufficiently slow translocation in the short nanopore, 
i.e., under low driving force and/or high friction, the ex- 
ponent a in the scaling law r A^" for the transloca- 
tion time T can be estimated by the balance of driving 
force F and frictional force. The velocity of the center of 
mass along the direction of the driving force is w ^ F/C 
The chain is not severely deformed during slow translo- 
cation processes, and the chain moves a distance of Rg 
during the translocation. Thus, the translocation time 
becomes t w Rg/v - N^+^^/F 0, [13, H^- Similarly, 
for polymer translocation under a constant pulling force 
acting on one chain end, the polymer travels a distance 
of ~ N during translocation, and r « R\\/v ^ N^^/F 
[3^. However, if the external driving force is a function of 
time as for the present problem in a long channel with an 
intra-channcl force, the translocation dynamics becomes 
different, as outlined in the following. 

For a polymer of chain length N in a good bulk sol- 
vent in two dimensions (2D), the radius of gyration of 
the chain Rg scales as Rg ^ N'^^^a, where 1^20 ~ 0.75 is 
the Flory exponent in 2D, and a is the segment length. 
For a polymer confined between two strips embedded 
in 2D, the chain will extend along the channel to form 
blobs of size R, as long as i? > ct. Each blob contains 
g ~ (i?/cr)^/'^^" beads, and the number of blobs is Ub = 
N/g - nIu/RY/"^'' . The free energy cost for the chain 
confinement in units of ksT is T = N{a/Ry^''^'' [H]. 



The blob picture then predicts the longitudinal size of the 
chain to be = mR ^ Na{a/ Ry^"^"-^ ~ NR-^'^ 
[2^. The longitudinal relaxation time r|| is defined as 
the time needed by a polymer to move a distance of 
the order of its longitudinal size, . Thus, ry scales as 

T|| - - N^R-'^/^, where I) - 1/iV is the diffusion 

constant of the chain. 

For a completely confined polymer in a 2D long channel 
under an external field, the total driving force is NF with 
F being the force acting on one bead. If the chain moves a 
distance of the order of under weak driving forces, the 
time cost can be estimated as Tss ^ R\\/v ^ N R^^/^ / F , 
where v ^ F \s the velocity of the polymer. Under 
strong driving forces, R\\ ~ TV is still correct [45|, but 
^ R~^/^ breaks down due to the deformation of the 
blob and the elongation of the chain under the driving 
force. Thus, one can only obtain Tss ~ N/F. For poly- 
mer translocation through a long channel with length 
L ^ i?g, we therefore obtain t^s ^ L/v ^ 1/F, which is 
independent of TV [46| . This steady state behavior in the 
channel, valid sufficiently far away from the two channel 
ends, is due to the fact that the force F acts uniformly on 
each monomer. The behavior is different when wc con- 
sider the translocation time r measuring the time from 
threading the first monomer into the channel until the 
moment when the entire chain is inside the channel. 

To assess the capture of a polymer into the long chan- 
nel under the influence of the intra-channel field, the sit- 
uation is somewhat complicated due to the fact that the 
driving forces are proportional to the number of translo- 
cated beads. However, we still can give a rough estimate, 
if the chain monomers sufficiently quickly reach a steady 
state in the channel. If s{t) monomers are located inside 
the channel, the chain on the cis side experiences the fric- 
tion force Fcisj and the entropic force Fds^e, while the 
translocating part experiences the friction force Ftransj, 
the entropic force Ftrans,e and the driving force s(t)F. 
According to the balance of driving force, frictional forces 
and entropic forces, we obtain the force relation 

Fcis.f ~^ Fq^^.c ~I~ Fij^f^jig j^ -i~ -ftrans.e — s(^t^F. (4) 

Based on the blob picture, 

Ftransj « C^W^ « C R'^/' is{t)^ , (5) 

with C being a constant. For a chain with size TV = 128 
in a channel with R = 4.5 in units of a, the longitudinal 
size is i?|| K, 90.73, and then wc obtain C = 1.17. Once 
more than g monomers enter into the channel, Ftrans,e ~ 
2.12kBT/R according to field-theoretical methods for 2D 
geometries (strips) and hard- wall boundaries Un- 
der the assumption that the chain on the cis side is 
close to its equilibrium state and is not severely deformed 
{Fcis,eRg{N — s{t)) < ksT), Fds^c&n be estimated as 
F„s,e « kBTL^/[Rl{N - s{t))] 1^, with being the 
elongation of the chain along the x direction. In ad- 
dition, Fcis.f = £,^iLs{t)+i'^i^ which under sufficiently 
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slow translocation is negligible compared with Ftransj, 
due to the much slower velocity. Thus, at the beginning 
of the translocation, Fcis,e plays an important role. How- 
ever, once more than g monomers enter the nanopore, 
Fcis,e can be neglected compared with Ftrans,e- Due 
to the fact that the driving force s{t)F increases over 
time, after reaching a critical s{t) the translocation dy- 
namics is dominated by the balance between the driving 
force s{t)F and Ftransj, particularly for larger F where 
Ftrans.e bccomcs negligible. Then the translocation dy- 
namics is dominated by 



ds{t) 
dt 



s{t)F, 



such that 



(6) 



(7) 



If Ftrans,e cannot be neglected, once more than g 
monomers enter into the channel, the translocation dy- 
namics is controlled by 



CR 



-1/3 



ds{t) 
dt 



■Ft, 



s{t)F, 



(8) 



such that 



CR 



-1/3 



N CR-^'^iFu 



F 



F2 



'-[l~\n{FN -Ftrans.e)]- 

(9) 

The results in Eqs. ^ and ([9|) are based on Eq. ([5|), 
where the equilibrium value of is used. As mentioned 
above, under strong driving forces ~ N is still cor- 
rect, but i?|| ~ R~^/^ breaks down. Therefore, we do not 
expect the scaling exponents of t with R predicted in 
Eqs. ([7]) and © to be strictly correct. However, we can 
still use Eq. © with Ftrans,e ~ 2.12kBT/R for the equi- 
librium state to estimate the scaling of r with N. Fig. [2] 
demonstrates that the scaling t N for long chains is 
fulfilled, and Ftrans,e is thus important only for short 
chains. 

If during the the translocation process the monomers 
in the channel cannot sufficiently quickly assume steady 
state motion, the translocation dynamics is more compli- 
cated, as discussed in the following based on numerical 
results. 



B. simulation results 

We now proceed to present our numerical results and 
discuss them in the light of above scaling arguments. 

1. Translocation probability as function of the driving force 
F and channel width R 

In our simulations initially the first monomer is held 
in place right at the channel entrance. After equilibra- 
tion of the chain with fixed first monomer, the constraint 
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FIG. 2: Translocation time r predicted in Eq. ((9| is plotted 
as function of chain length A'^ for different channel widths R 
and driving forces F . 
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FIG. 3: Translocation probability as function of the driving 
force F for fixed channel width R = 4.5. The results for vari- 
ous chain lengths approximately coincide. The line indicates 
a linear relation between the translocation probability and 
the driving force F per beads. 



is relieved. At high driving force F ^ already acting on 
the first monomer, the probability is relatively high that 
the chain will actually completely move into the chan- 
nel, and not escape back to the cis side. At lower F 
retraction becomes more frequent. In this context we re- 
fer to a successful translocation as the event when the 
chain fully enters the pore, i.e., all N monomers are in- 
side the channel such that the chain reaches the translo- 
cation coordinate s = N. If the first monomer retracts 
from the channel, the translocation is viewed unsuccess- 
ful. In Fig. 13] we show the translocation probability 
of successful events as function of the driving force F, 
and for different chain lengths N. For each value of F 
and N we perform as many translocation attempts until 
1000 successful translocation events are reached. In the 
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FIG. 4; Translocation probability as function of the channel 
width R for driving forces F — 0.2 and 0.5, and for two chain 
lengths. 
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for different channel widths R and driving forces F. 



investigated force range the translocation probability in- 
creases almost linearly with increasing translocation force 
F, and is approximately independent of the chain length 
TV. This implies that it is important to capture the first 
few monomers in the channel, such that the overall force 
s{t)F reaches a sufficiently large value to prevent chain 
retraction from the pore. 

Fig. m depicts the translocation probability as func- 
tion of the channel width R for two chain lengths, at 
driving forces F = 0.2 and 0.5. The resulting curves 
show an increase from quite small channel width towards 
wider channels, as expected. This is due to the decreased 
entropic resistance against threading of the chain into 
the channel. The approximate independence of the chain 
length suggests that this is, again, a more local effect: it 
matters to succeed sucking the first few monomers into 
the channel to ensure complete translocation. At larger 
values for R we observe a flattening of the curve, likely 
due to the elongation of the chain in the channel. 

2. Translocation times as function of the chain length N 
and channel width R 

Fig. [5]shows the translocation time r as function of the 
chain length N for different driving forces F and channel 
widths R. For > 32, a good scaling behavior r ^ N" is 
observed, with scaling exponent a « 1 for F = 0.1 ... 0.5 
and R = 4.5, as well as for F = 0.5 and R ~ 6.0. In 
these cases, the driving force dominates the transloca- 
tion dynamics, and the value of a is in agreement with 
the result ([T]) obtained from our scaling arguments. For 
N < 32, Ftrans.e ~ 2.\2kBT / R incrcascs and is no longer 
negligible for weak driving forces as predicted in Eq. Q 
and plotted in Fig. [5J Decreasing R for fixed F = 0.2, 
the scaling exponent a is reduced to 0.88 ± 0.01, and 
0.81 ± 0.01 for R = 3.0 and 2.0, respectively. The relax- 
ation time of the chain increases with decreasing i?, which 



leads to significant non-equilibrium situations. This fac- 
tor leads to a more complicated translocation dynamics, 
that can no longer be grasped by our simple scaling argu- 
ments leading to Eqs. ([7]) and ©. The resulting behav- 
ior, however, is good news: increasing the overall force 
exerted on the chain during the process of threading the 
chain into the channel leads to a more efficient transloca- 
tion, with a scaling exponent smaller than one. It should 
be noted that once the entire chain is fully in the chan- 
nel and reaches a steady state, the typical time to cover a 
distance R\\ becomes independent of iV, as shown above. 

The radius of gyration of the polymer along the chan- 
nel axis X, which is proportional to , at the moment 
just after complete entrance into the channel, for dif- 
ferent values of the translocation force F and channel 
width R turns out to follow the almost linear relation 
Rg.x ^ N. This is demonstrated in Fig. [HI This is 
perfectly in line with our scaling arguments for devel- 
oped above. The inset in Fig. [B] shows the translocation 
velocity, which decreases with chain length N, and al- 
most saturates for larger N. With increasing channel 
width i?, the translocation time decreases rapidly and 
then saturates for larger R, as shown in Fig. [T] Likely 
this behavior is due to the decaying influence of the en- 
tropic force in the channel, Ftrans,e ^ kBT/R and the 
elongation of the chain under the driving force. 



S. Translocation time as function of the driving force F 

Fig. [S]shows the translocation time r as function of the 
driving force F for different chain lengths A^. We observe 
a behavior close to the predicted inverse proportionality 
with F, T ^ F^^ from Eq. ([5]), for different chain lengths 
when the channel entropic force Ftrans,e ^ ksT / R is less 
pronounced. However, for smaller R, such as i? = 2.0, 
the entropic force Ftrans,e plays a more important role 
in the translocation dynamics: thus, for TV = 64 and 
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FIG. 6: Radius of gyration of the translocating polymer at 
the moment just after translocation, i.e., full entrance into 
the channel, for different translocation forces F and channel 
widths R. The insert shows the translocation velocity d as a 
function of the chain length A^. 



1200 
1000 

800 
600 

400 



200 



1 1 1 1 1 1 1 1 




-□-N=128, F=0.2 1 






-o-N=128, F=0.5 ■: 






~~Q 5 ; 












§ Q \ 


1 1 1 1 1 1 1 1 
2 


3 


1 1 1 

R 4 5 6 7 




100 150 

s 

FIG. 9: Waiting time distribution for different chain lengths 
under driving force F = 0.2 and with channel width R — 4.5. 
We define the waiting time of bead s as the average time 
between the events that bead s and bead s + 1 move away 
from the channel entrance. 
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FIG. 7: Translocation time r as function of the channel width 
i? for iV = 128, and F = 0.2 and 0.5. 



FIG. 10: Waiting time distribution under different driving 
forces for chain length A'^ = 128 and channel width R = 4.5. 
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FIG. 8: Translocation time r as function of the driving force 
F. 



R = 2.0, the chain cannot fully enter the nanopore at all, 
while for F > 0.2. we observe r F~^-^^. 



4. Waiting time distribution 

The dynamics of a single segment passing through the 
pore during translocation is an important quantity con- 
siderably affected by different driving mechanisms. The 
nonequilibrium nature of the translocation process has 
significant effects on this quantity. We numerically cal- 
culated the waiting times for all beads in a chain of length 
A^. The waiting time of bead s for successful transloca- 
tion hereby is defined as the average time between the 
events that bead s and bead s + 1 move away from the 
channel entrance and further into the channel. 

Fig. [9] shows the waiting time distribution for different 
chain lengths under driving force = 2.0 and with chan- 
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FIG. 11: Waiting time distribution for different channel width 
R and chain length N = 64 under driving force F = 0.2. 



nel embedded in a two dimensional geometry, mimick- 
ing typical nano/micro-fluidic setups used, for instance, 
for DNA analysis. We analyzed how the translocation 
dynamics depends on the chain length N of the chain- 
to-be-translocated, the driving force F, and the channel 
width R. During the the translocation process if the 
monomers in the channel can sufficiently quickly assume 
steady state motion, we observe the scaling t ^ N/F 
of the translocation time r with F and TV. With smaller 
channel width R, steady state motion cannot be achieved, 
effecting a non-universal dependence of r on A'^ and F. 
We also find that the waiting time distribution shows a 
maximum at the same translocation coordinate Smax for 
different chain lengths, and Smax decreases with increas- 
ing driving force or channel width. We believe that this 
study opens up a new aspect to translocation studies, 
adding relevant quantitative information on the growing 
field of biopolymer analysis in nanochanncls. 



nel width R = 4.5. Interestingly, the maximum waiting 
time occurs at Smax ~ 10 for all measured chain lengths. 
With increasing s, the waiting time increases rapidly to 
the maximum value, before a significantly softer decay. 
For increasing chain length a plateau appears. This form 
of the waiting time distribution again suggests that the 
initial capture of the chain in the channel is a local pro- 
cess, unaffected by the remainder of the chain on the cis 
side. At later stages of the process the leftover monomers 
are sucked into the channel more easily, due to the high 
force s{t)F reached by the s(t) monomers in the chan- 
nel. Fig. [TU] shows the waiting time distribution un- 
der different driving forces F, for chain length A^ = 128 
and channel width R = 4.5. We observe that the num- 
ber of monomers corresponding to the maximum waiting 
time, Smax, shifts to smaller values with increasing driv- 
ing force, as one would expect. Increasing R, we find 
that Smax rapidly decreases but saturates for larger R, 
compare Fig. [TTl 

The coming into existence of non-equilibrium situa- 
tions in the translocation process, we display typical 
chain configurations in the Appendix. Due to the in- 
creasing driving force s{t)F, the chain on the cis side 
experiences trumpet, stem-flower, and straight chain con- 
formations during the translocation process, see Figs. [121 ^ 

M 

IV. CONCLUSIONS 

Using Langcvin dynamics simulations, we investigated 
the dynamics of polymer translocation into a long chan- 
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Appendix: Translocation snapshots 

To get an idea of translocation process, we show the 
typical chain conformation for different stages of the 
translocation process. The chain length A^ = 128. Figs. 
1121 and [T3] shows chain conformation during transloca- 
tion for channel width R = 4.5 under the driving force 
F — 0.1 and F = 0.5, respectively. For F = 0.1, 
at t = 100 the typical trumpet chain conformation oc- 
curs. At t = 100 and t = 1200, we observe the typ- 
ical stem-flower conformation. Finally, the chain is al- 
most straight a.t t = 1500. For F = 0.5, the chain is 
stretched to straight conformation at t = 300. Figs. [HI 
shows chain conformation during translocation for chan- 
nel width i? = 2.0 under the driving force F = 0.2. The 
chain in the channel is very straight during the translo- 
cation. 
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